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ABSTRACT Some features of hydro- and thermodynamics, as applied to atmospheres 
and to stellar structures, are puzzling: 1. The suggestion, first made by Laplace, that our 
atmosphere has an adiabatic temperature distribution, is confirmed for the lower layers, 
but the reason why it should be so is understood only qualitatively. 2. Arguments in which 
a concept of energy plays a role, in the context of hydro-thermo dynamical systems and 
gravitation, are often flawed, and some familiar results concerning the stability of model 
stellar structures, first advanced at the end of the 19th century and repeated in the most 
modern textbooks, are less than completely convincing. 3. The standard treatment of 
relativistic thermodynamics does not allow for a systematic treatment of mixtures, such 
as the mixture of a perfect gas with radiation. 4. The concept of mass in applications of 
general relativity to stellar structure is unsatisfactory. It is proposed that a formulation 
of thermodynamics as an action principle may be a suitable approach to adopt for a new 
investigation of these matters. 

We formulate thermodynamics of ideal gases in terms of an action principle and study 
the interaction between an ideal gas and the photon gas, or heat. The action principle 
provides a hamiltonian functional, not available in traditional approaches where familiar 
expressions for the energy have no operative meaning. The usual polytropic atmosphere 
in an external gravitational field is examined, in order to determine to what extent it 
is shaped by radiation. It is easy to understand that radiation sustains the atmosphere 
(prevents cooling) , but the temperature profile is largely determined by intrinsic properties 
of the gas and it is difficult to interpret it as an effect of radiation. This has led some 
people to question whether the rule of uniform temperature as an absolute condition for 
equilibrium is valid in the presence of a gravitational field. An experiment that involves a 
centrifuge and that has wider implications in view of the equivalence principle, is proposed, 
to ascertain the influence of gravitation on the equilibrium distribution with a very high 
degree of precision. A new formulation of the concept of radiative equilibrium is proposed. 

The choice of boundary conditions for radial, stellar stability calculations is clarified 
with the help of a properly defined, conserved mass distribution. 

PACS Keywords: Atmosphere, photon gas, action principle. 
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I. Introduction 

The statement that any two thermodynamic systems, each in a state of equilibrium 
with a well defined temperature, in thermal equilibrium with each other, must have the 
same temperature, is a central tenet of thermodynamics. A natural extension is that the 
temperature in an extended but closed system in a state of equilibrium must be uniform, 
but there does not seem to be universal agreement on whether this is true in the presence 
of gravitational fields. The question comes up in the investigation of terrestrial or stellar 
atmospheres, where the gravitational forces create a non-uniform density distribution. 

This paper is a study of the polytropic atmosphere. We want to know if (or to what 
extent) the polytropic relations are to be attributed to intrinsic properties of the gas, or 
to radiation. In this we are guided by a strong preference for action principles. 

The textbooks present hydrodynamics as the theory of a continuous distribution of 
matter, described in the simplest case by two fields or distributions: a density field and a 
velocity field, both defined over 1R 3 or a portion thereof. The role of temperature is often 
disguised, assumed to be determined by the density and the pressure. Classical thermo- 
dynamics, on the other hand, is primarily the study of states of equilibrium, with uniform 
density and temperature, and relations between such states. In this context extremum 
principles play a pivotal role; see for example Callen (1960). Texts that deal with flow 
of matter or with temperatures that vary in time or from one point to another are found 
under the heading of heat transfer, fluctuations or thermodynamics of irreversible pro- 
cesses. See for example Stanyukovich (1960), Miiller (2007). These studies rely heavily on 
conservation laws, but variational principles are rarely mentioned. Texts that most closely 
parallell the present work are found under radiation hydrodynamics (Castor 2004). 

In this introduction we study a simple system from the point of view of hydrodynamics, 
on the basis of a well known action principle. The concept of temperature appears, but not 
as a dynamical variable. We offer a brief review of the history of the polytropic atmosphere 
(Section 1.6) and stress the role of mass. In Section II, we extend the action principle to 
include the temperature as an independent field variable. A lagrangian describes the 
configurations on a single polytrope. The Euler-Lagrange equations include the gas law as 
well as the polytropic relation between pressure and density. The familiar expression for 
the internal energy of an ideal, polytropic gas coincides with the hamiltonian; this appears 
to be a significant result. The energy and the pressure of radiation are incorporated in a 
natural way (Section II. 6). 

The theory is mathematically complete in the sense that no additional input from 
underlying microscopic physics is needed; as an example we derive a virial theorem that is 
proper to the action principle (Section II. 7). But it is physically incomplete since the role 
of convection and heat flow is not explicitly accounted for. 

To contribute to the debate on the question of whether an isolated atmosphere in a 
gravitational field tends to isothermal equilibrium we study an ideal gas in a centrifugue 
and invoke the equivalence principle (Section II. 8). 

In Section III we study the effect of sources of heat that induce transitions between 
polytropes. The radiation field contributes as a source of entropy. The principal idea be- 
hind this work was the hope of understanding, in a quantitative way, why real atmospheres 
tend to be polytropic; in this we have had only limited success. 
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In Section IV we take up the problem of the stability of polytropic, atmospheric 
models. A principal advantage of the method is that it provides us with a hamiltonian, 
expressed in terms of the dynamical variables. Some classical stability studies are found 
wanting, because of ad hoc definitions of various energies, and inapropriate boundary 
conditions. 

Section V makes the passage to General Relativity. Section VI has a summary of 
conclusions and several proposals for additional work, theoretical as well as experimental. 



1.1. Hydrodynamics 

The textbook introduction to hydrodynamics deals with a density field p and a velocity 
field v over R 3 , subject to two fundamental equations, the equation of continuity, 

p + div(H=o, p-=% (i-i) 

and the hydrodynamical equation (Bernoulli 1738) 

— grad p = p^-v := p(v + v • grad v). (1.2) 

This involves another field, the scalar field p, interpreted as the local pressure. The theory 
is incomplete and requires an additional equation relating p to p. It is always assumed 
that this relation is local, giving p(x) in terms of the density at the same point x, and 
instantaneous. 



1.2. Laminar flow 

Since we are reluctant to take on difficult problems of turbulence, we shall assume, 
here and throughout, that the velocity field can be represented as the gradient of a scalar 
field, 

v = —grad (1.3) 

In this case the hydrodynamical condition is reduced to 

grad p = p grad ($ - v 2 /2). (1.4) 

To complete this system one needs a relation between the fields p and p. 
Assume that there is a local functional V[p] such that 

p = p V' - V, V := dV/dp. (1.5) 

In this case dp = p dV and the equation becomes, if p ^ 0, 

grad V = grad ($ - v 2 /2) (1.6) 

or 

V' = $ - v 2 /2 + A, A constant. (1.7) 

The potential V[p] is defined by p modulo a linear term, so that the appearance of an 
arbitrary constant is natural. It will serve as a Lagrange multiplier. 

The introduction of a velocity potential guarantees the existence of a first integral of 
the motion, a conserved energy functional that will play an important role in the theory. 



3 



1.3. Variational formulation 

Having restricted our scope, to account for laminar flows only, we have reduced the 
fundamental equations of simple hydrodynamics to the following two equations, 



<9 M J M = 0, J* := p, J := pv 
V' = $ - v 2 /2 + A, 



, 2 . . (1-8) 



together with the defining equations 

v = -grad p := pV - V. (1.9) 

It is well known that these equations are the Euler-Lagrange equations associated with the 
action (Fetter and Walecka 1980) 

A = Jdtd 3 xC, £ = p(<5> -v 2 /2 + \) -V[p\. (1.10) 

The value of this last circumstance lies in the fact that the variational principle is a much 
better starting point for generalizations, including the incorporation of symmetries, of spe- 
cial relativity, and the inclusion of electromagnetic and gravitational interactions. It also 
gives us a valid concept of a total energy functional. 

1.4. On shell relations 

The action (1.10) contains only the fields $ and p, and the Euler-Lagrange equations 
define a complete dynamical framework, but only after specification of the functional V[p]. 
The pressure was defined by Eq.(1.9), p := pV — V, and one easily verifies that, on the 
trajectory, by virtue of the equations of motion, 

p = C (on shell). (1.11) 

This fact has been noted, and has led to the suggestion that the action principle reduce to 
the minimization of J p with respect to variations of p defined by thermodynamics (Taub 
1954), (Bardeen 1970), (Schutz 1970). But more is needed, including an off shell action. 
After adopting the action (1.10) it remains to relate the choice of the potential V to the 
thermodynamical properties of the fluid. 



1.5. Equation of state and equation of change 

An ideal gas at equilibrium, with constant temperature, obeys the gas law 

p/p = TIT. (1.12) 
Pressure and density are in cgs units and 

K = (1/aO x .8314 x 10 8 erg/K, 
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where fi is the atomic weight. In this paper the validity of the gas law is assumed to 
hold, locally at each point of the gas, under all circumstances, including the case that 
gravitational and electromagnetic fields are present. Effective values of fx are 

Atomic hydrogen : p = 1, Air : p = 29, Sun : fx = 2. 

Equation (1.12) is the only equation that is referred to as an 'equation of state'. Other 
relations, to be introduced next, are 'equations of change', this term taken from Emden's 
"Zustandsanderung" , for their meaning is of an entirely different sort. The most important 
is the polytropic relation 

p = Ap 1 \ A, 7' constant. (1-13) 

This relation defines a polytropic path or polytrope in the p, v diagram (v = 1/p). A 
polytropic atmosphere is one in which, as one moves through the gas, the variables p and 
p change so as to remain always on the same polytrope, the temperature being determined 
by Eq.(1.12) always. Eq.(1.13) is a statement about the system, not about the gas per se. 
The index of the polytrope is the positive number n defined by 

n 

Important special cases are 

n' = 0, 7 = oo, p = constant, 

y = Cp/Cv, specific entropy = constant, 
n' = oo, 7' = 1, T = constant. 
The numbers 7, n are defined by 

1 



7 : = Cp/C v =: 1 + 



n 



The number n is the adiabatic index of the gas. According to statistal mechanics 2n is the 
number of degrees of freedom of each molecule in the gas. That atmospheres tend to be 
polytropic is an empirical fact. 

The case that 7' = 7 is of a special significance. A polytrope with 7' = 7 is a path 
of constant specific entropy; changes along such polytropes are reversible and adiabatic; 
these polytropes and no others are adiabats. 

Fix the constants A, 7' and consider an associated stationary, polytropic atmosphere. 
Since both (1.12) and (1.13) hold we have 

p = const, p 1 , p = const. TV- 1 , T = const, p 1-1 ^ 1 . (1-14) 

In any displacement along a polytrope from a point with pressure p and temperature T, 
we shall have dp/ p = (l/Y)dp/p, so that an increase in pressure leads to an increase in 
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density that is greater for a smaller value of 7'. If a parcel of gas in this atmosphere is 
pushed down to a region of higher pressure, by a reversible process, then it will adjust 
to the ambient pressure. If 7 > 7', then it will acquire a density that is lower than the 
environment; it will then rise back up; this atmosphere is stable. But if 7' > 7 then the 
parcel will be denser than the environment and it will sink further; this atmosphere is 
unstable. Thus we have: 

A stable, polytropic atmosphere must have 7' < 7, n' > n. 

Most stable is the isothermal atmosphere, 7' = 1. 

In hydrodynamics, the isothermal atmosphere can be given a lagrangian treatment by 
taking 

V = KTp\ogp. (1.15) 

We suppose that the gas is confined to the section zq < z < zo + h of a vertical cylinder 
with base area A and expect the density to fall off at higher altitudes. A plausible action 
density, for a perfect gas at constant temperature T in a constant gravitational field 4> = gz, 
g constant, is 

£[$, p ] = p ($ - v 2 /2 -gz + X)- TLTp log p. (1.16) 

We may consider this an isolated system with fixed mass and fixed extension. 

At equilibrium $ = 0,t7=0,p = and the equation of motion is V = TIT (1 + log p) = 
X — gz, hence 

P {X, y, Z) = e -l+V*T e -gz/*T M = ^^I e -l+A/^T (1 _ e -gh/R T) e -gz„/^T 

g 

and after elimination of A 

gM e -g( z - z o)/^T gjy[ e -g(z-z )/TCT 

9 = ART 1 - e-sVRT ' p = ~A 1 _ e -gh/RT • ( L17 ) 

There is no difficulty in taking the limit h — > 00. The volume becomes infinite but it can be 
replaced as a variable by the parameter zq. This atmosphere is stable; a proof is presented 
in Section IV. 1. 

The isothermal atmosphere is usually abandoned in favor of the polytropic atmo- 
sphere. 

A polytropic gas can be described by the lagrangian (1.10), with 

V = ap 1 , a, 7' constant. 
Variation with respect to p gives 

a 7' 1 '1 

u u 
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to be supplemented by the gas law, Eq.(L12). Among the many applications the following 
are perhaps the most important. In the case of sound propagation the gas is initially 
awakened from equilibrial turpor and then left in an isolated, frenzied state of oscillating 
density and pressure, with the temperature keeping pace in obedience to the gas law. All 
three of the relations (1.14) are believed to hold, with 7' = 7. The oscillations are usually 
too rapid for the heat to disseminate and equalize the temperature, so that the neglect of 
heat transfer may be justified. In applications to the atmospheres one uses the polytropic 
equation of change (1.13) and obtains the temperature from the gas law. Understanding 
the resultant temperature gradient in terms of convection, or as the effect of the heating 
of the air by solar radiation, or both, is one of the main issues on which we had hoped to 
gain some understanding. 

At mechanical equilibrium v = 0, p = and A — gz = a^p 1 / n , hence 

a7 

Since the density must be positive one does not fix the volume but assumes that the 
atmosphere ends at the point z\ = X/g. Then 

M = A(^-) n f Z \z 1 - z) n dz = _J^_(^)n. 
a 7 ' J Z0 n + l v a 7 ; 

This fixes h and thus z\ and A. If the atmosphere is an ideal gas then the temperature 
varies with altitude according to 

nT = P / P = V /n = s^ (1.18) 

n ri+1 

Because the lagrangian does not contain T as a dynamical variable it is possible to impose 
this condition by hand. 

One would not apply this theory down to the absolute zero of temperature, but even 
without going to such extremes it seems risky to be predicting the temperature of the 
atmosphere without having made any explicit assumptions about the absorption or gener- 
ation of heat that is said to be required to sustain it. Yet this has been the basis for the 
phenomenology of stellar structure, as well as the earth's atmosphere, from the beginning 
(e.g. Lane 1870, Ritter 1878). 

The success of the isotropic model is notorious, and this success can be explained in 
physical terms, but the theory is incomplete since it does not account for heat flow, nor 
convection, both of which are needed to compllete the picture. 

For air, with atomic weight 29, 1Z = 2.87 x 10 6 ergs/gK and n = 2.5. At sea level, 
g = 980cm/sec 2 , the density is p = 1.2 x 10~ 3 g/cm 3 , the pressure p = 1.013 x 10 6 dyn/cm 2 . 
Thus 

p/p = .844 x 10 9 cm 2 /sec 2 , T = T = 294K, z x = 3.014 x 10 6 cm « 30km. 

and the dry lapse rate at low altitudes is — T' = 294/ z\ = 9.75K/km. The opacity that 
is implied by this is mainly due to the presence of CO2 in the atmosphere. Humidity 
increases the opacity and decreases the lapse rate by as much as a factor of 2. 



7 



1.6. Historical notes on polytropic atmospheres 

Observations of reversible transformations of near-ideal gases, carried out during the 
19th century, can be summarized in what is sometimes called the laws of Poisson, 

' ' ' 1 

p(xT n \ pocT n ' +1 , pocp 1 ', 7 ; = H constant. 

n' 

In the original context all the variables are constant and uniform. The exponents as 
well as the coefficients of proportionality are the same for all states that are related by 
reversible transformations. Statistical mechanics explained this result and confirmed the 
experimental value 7' = 7 = Cp/Cy. As far as can be ascertained, the presence of 
terrestrial gravitation and ambient radiation had no effect on these experiments. In a first 
extrapolation the same relations were taken to hold locally in dynamic situations, as in 
the case of sound propagation. The gas is not in equilibrium and the variation of the 
temperature from point to point, and with time, is obtained from the gas law. 

For the atmosphere of the earth it was at first proposed that the temperature would 
be uniform. However, the existence of a temperature gradient was soon accepted as an 
incontrovertible experimental fact. The first recorded recognition of this, together with an 
attempt at explaining the same, may be that of Carnot, in the paper in which he created 
the science of thermodynamics (Carnot 1824). Carnot quotes Laplace: "N'est-ce pas au 
refroidissement de l'air par la dilatation qu'il faut attribuer le froid des regions superieures 
de l'atmosphere? Les raisons donnees jusqu'ici pour expliquer ce froid sont tout a fait 
insuffisantes; on dit que Pair des regions elevees, recevant peu de chaleur reflechie par la 
terre, et rayonnant lui meme vers les espaces celestes, devait perdre de calorique, et que 
c'etait la la cause de son refroidissement; ... " This may be the first time that the influence 
of radiation is invoked. The temperature gradient is attributed to the greenhouse effect, 
and Laplace was an early skeptic, for he continues "...mais cette explication ce trouve 
detruite si Ton remarque qu'a egale hauteur le froid regne aussi bien et meme avec plus 
d'intensite sur les plaines elevees que sur les sommets des montagnes ou que dans les parties 
d'atmosphere eloignees du sol." It is not clear that the two explanations are at odds with 
each other; Laplace apparently postulates that the atmospheres over lands at different 
elevations are related by adiabatic transformations, without explaining why. 

By rejecting the role of radiation as the cause of the temperature gradient, Laplace 
seems to suggest that the same would be observed in an atmosphere subject to gravitation 
but totally isolated from radiation, neither exposed to the radiation coming from the sun 
nor radiating outwards. As was strongly emphasized in the later phases of this debate, this 
would contradict the belief that the thermal equilibrium of any isolated system, gravitation 
and other external forces notwithstanding, is characterized by a uniform temperature. 

In 1862 W. Thomson, in the paper "On the convective equilibrium of the temperature 
in the atmosphere" , defines convective equilibrium with these words "When all parts of a 
fluid are freely interchanged and not sensibly influenced by radiation and conduction, the 
temperature is said to be in a state of convective equilibrium." He then goes on to say that 
an atmosphere that is in convective equilibrium is a polytrope, and we think that he means 
an adiabat, since this is probably implied by the words "freely interchanged", although 
the value of 7 is taken from experiment and not from statistical mechanics. At first sight 
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the clause "and not sensibly influenced by radiation" would seem to imply that an isolated 
atmosphere has a temperature gradient, but this conclusion would be premature, as we 
shall see. 

In 1870 H.J. Lane made the bold assertion that the laws of Poisson may be satisfied 
in the Sun. The terrestrial atmosphere (or part of it) had already been found to be 
well represented by the same relations. Referring to Lane's paper Thomson, now lord 
Kelvin, explains how convective equilibrium comes about (Thomson 1907). He argues 
that the atmosphere is not, cannot be, at rest, and this time radiation plays an essential 
role. The upper layers loose heat by radiation and the lower temperature leads to an 
increase in density. This produces a downward current that mixes with a compensating 
upward drift of warmer air. This continuing mixing takes place on a time scale that is 
too short for adjacent currents to exchange a significant amount of heat by conduction 
or radiation, especially since the variations of temperature are very small. It is evident 
that Thomson offers his explanation of the temperature gradient to account for its absence 
in an isolated atmosphere, for he says that, "an ideal atmosphere, perfectly isolated from 
absorption as well as emission of radiation, will, after enough time has passed, reach a state 
of uniform temperature, irrespective of the presence of the gravitational field" . Thomson 
accepts the mechanism of Laplace and Carnot, as it is at work in the real atmosphere, 
but he goes further. He believes that the lower temperature aloft is intimately tied to the 
existence of radiation, implying that it is driven by net outwards radiation. (The effect 
of solar radiation on the terrestrial atmosphere is not explicitly mentioned.) It is difficult 
to tell whether or not Thomson is in disagreement with Laplace, but the precision of his 
statements represents a marked improvement over his predecessors and his earlier work. 

The principal developers of the field, Ritter (1878-1883) and Emden (1907), seem 
to accept the idea of convective equilibrium. It may be pointed out, however, that this 
mechanism is in no way expressed by the equations that these and other authors use 
to predict the behaviour of real atmospheres. The concept of convective equilibrium is 
introduced to one purpose only: to avoid contradiction with the ideas on thermal equilibrium 
of isolated systems. It receives no quantitative theoretical treatment. 

Nor was it accepted by everybody. A famous incidence involves Loschmidt (1876), 
who believed that an isolated atmosphere, at equilibrium in a gravitational field, would 
have a temperature gradient. But arguments presented by Maxwell and Boltzmann (1896) 
led Loschmidt to withdraw his objections, which is hardly surprising given the authority 
of these two. Nevertheless, it may be pointed out that no attempt was made, to our 
knowledge, to settle the question experimentally. See however (Graeff 2008). 

An alternative to convective equilibrium was proposed by Schwarzschild (1906) and 
critically examined by Emden. To understand how it works we turn to Emden's book 
of 1907, beginning on page 320. Here he invokes a concept that is conspicuously absent 
from all his calculations on polytropic spheres in the rest of the book: heat flow. It must 
be agreed that the atmosphere is not completely transparent, and that heat flow is an 
inevitable consequence of the existence of a temperature gradient. The most important 
observation is that heat flow is possible in stationary configurations (T = 0) provided that 
the temperature gradient is constant. The heat flux due to conduction and radiation can 
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be expressed as 

F = Cx/T, F i = C ij d j T, 

where the tensor C includes the thermal conductivity as well as the effective coefficient 
of heat transfer by radiation. The divergence of the flux is the time rate of change of 
the temperature due to conduction and radiation. In a stationary, terrestrial atmosphere, 
with no local energy creation, this must vanish. Emden's atmospheres are polytropes, 
with temperature gradients that are constant (it appears that he takes C to be constant) . 
That is interesting, for it reminds us that the entire edifice implicitly demands that this 
condition, of a constant temperature gradient, be satisfied. 

We note that the direction of flow is from hot to cold, outwards. Confining ourselves to 
planetary atmospheres, with no local energy generation, this calls for an explanation, since 
the ultimate source of energy is above. Here we have to return to the oldest explanation of 
the existence of a temperature gradient, dismissed by Laplace (op. cit.): the greenhouse 
effect. The atmosphere is highly transparent to the (high frequency) radiation from the 
Sun but opaque to the thermal radiation to which it is converted by the ground. The 
atmosphere is thus heated from below! 

If the atmosphere is stable in the sense discussed above, when 7' < Cp/Cy, then it is 
not necessary to assume that any convection takes place. In this case one speaks of (stable) 
'radiative equilibrium'. Convective equilibrium steps in when the stationary atmosphere 
is unstable, but it is no longer used to explain the existence of a temperature gradient; it 
is the effect rather than the cause of it. 

A difficulty is present in all accounts of stellar structure up to 1920. The energy ob- 
served to be emitted by the Sun, attributed to contraction of the mass and the concomitant 
release of internal energy, was far too small to account for the age of the sun as indicated 
by the geological record. The situation changed with the discovery of thermonuclear en- 
ergy generation. Now there is plenty of energy available. At the same time there arose the 
realization that convection sometimes plays a very modest role; the concept of convective 
equilibrium was put aside and with it, Kelvin's explanation of the temperature gradient. 
According to Eddington (1926), "convective equilibrium" must be replaced by "radiative 
equilibrium" . He does not claim that this new concept explains the temperature gradient 
as well as Kelvin's convective equilibrium does, but in fact the local generation of heat 
by thermonuclear processes creates an outward flow of heat and a negative temperature 
gradient. 

Emden's implicit invocation of the heat equation reminds us that this equation prob- 
ably should replace the simple rule that 'a system in equilibrium must have a uniform 
temperature'. 

Finally, there emerges a physical picture that seems to account for all the principal 
features of some atmospheres, in a qualitative way. If one begins with an isolated atmo- 
sphere in equilibrium one may take it as an axiom that the temperature is uniform. A 
weak dose of radiation upsets the static equilibrium and the gas goes into a state of con- 
vective equilibrium, characterized by an adiabatic temperature profile. When the intensity 
of the radiation increases it produces an outward heat flux and with it an added pressure. 
This pressure is the gradient of the energy of radiation and the additional 'potential' has 
to be overcome when a parcel of gas is moving downward; that is, the polytropic index 
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is increased. This in turn stabilizes the atmosphere against convection and eventually it 
reaches a state of 'radiative equilibrium', with insignificant convection. See (Cox and Giuli 
1968), page 271. If this interpretation is correct then the theory must be completed by 
inclusion of heat flow as an additional degree of freedom, and by an account of the prop- 
erties of the coefficient C. A more distant goal would be the inclusion of convection into 
the mathematical desription. 



1.7. The mass 

To speak of a definite, isolated physical system we must fix some attributes, and 
among such defining properties we shall include the mass. We insist on this as it shall turn 
out to be crucial to the logical coherence of the theory. The density p will be taken to 
have the interpretation of mass density, and the total mass is the constant of the motion 

M = J d 3 x p. 

Such integrals, with no limits indicated, are over the domain E of definition of p and is the 
total extension of our system in R 3 . 

Since the total mass is a constant of the motion it is natural to fix it in advance and 
to vary the action subject to the constraint J E d 3 x p(x) = M. We introduce a Lagrange 
multiplier and the action takes the form 

A = J d 3 x( K p(§-v 2 /2)-V^+\( y J d 3 xp-My (1.19) 

In the simplest case of a polytropic equation of state and no external forces we get the 
following equations of motion 

^ = 0, V' = $- v 2 /2 + A = a 7 p 1/n . 

Here A is to be chosen for each solution so as to satisfy the constraint. In the case of a 
static solution with $ = 0, v = the density is constant. Assuming a finite system with 
volume V we have M = pV = {\/d^) n V and since M is given, 

M a,M s „ , „ .M s1 /„ 

p=y> p= n i V ) ' A = ^(y) v ( L2 °) 

The conservation of mass has important implications for boundary conditions. 
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II. The first law 

II. 1. Thermodynamic equilibrium 

A state of thermodynamical equilibrium of a system that consists of a very large 
number of identical particles is denned by the values of 3 variables, a priori independent, 
the density D, the pressure P and the temperature T. These are variables taking real 
values; they apply to the system as a whole. In the case of any particular system there is 
one relation that holds for all equilibrium states, of the form 

T = f(D,P). 

It is written in this form, rather than F(T, D, P) = 0, because a unique value of T is 
needed to define a state of equilibrium between two systems that are in thermal contact 
with each other: it is necessary and sufficient that they have the same temperature. This 
statement incorporates the zeroth law. 

If we divide our system into subsystems then these will be in thermal equilibrium with 
each other only if they have the same temperature. This, at least, is the inherited wisdom; 
we shall honor it as long as possible. 

The ideal gas at equilibrium is defined by global variables T, D, P, and two relations. 
The principal one is the gas law 

P/D = TZT, TZ = .8314 x 10 8 ergs/K, 

where 1/D is the volume of a mole of gas. The other may take the form of an expression 
for the internal energy. 



II. 2. The ideal gas in statistical mechanics 

Here again we consider a gas that consists of identical particles (Boltzmann statistics), 
each with mass m and subject to no forces. It is assumed that the ith particle has momen- 
tum pi and kinetic energy p 2 /2m. This is an ideal gas, satisfying the relation P/D = TIT 
at equilibrium. It is assumed that the number N of particles with energy E is given by 
the Maxwell distribution 

N(E) oc e~ E / kT , (2.1) 

which implies a constant density in configuration space. Now place this gas in a constant 
gravitational field, with potential (p(x, y, z) = gz, g constant. Since the potential varies 
extremely slowly on the atomic scale it is plausible that, at equilibrium, each horizontal 
layer (</> constant) is characterized by a constant value of the temperature, density and 
pressure. Since neighbouring layers are in thermal contact with each other the temperature 
must (?) be the same throughout, 

T(z) = T = constant, 

and 

p(z)/p(z) = TIT, z> 0. (2.2) 
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The energy of a particle at level z is p 2 /2m + mgz and (2.1) now implies the following 
distribution in configuration space, 

p(x, y, z) oc e - m9Z ' kT , (2.3) 

in perfect agreement with (1.17). This supports the appearance of the logarithm in the 
expression for the potential, Eq.(1.15). Both derivations of the distribution rest on the 
assumption that the temperature is constant throughout the system. 

We conclude that the static solutions of the action principle, with action density 
(1.16) and T fixed, describe the equilibrium states of an ideal gas at fixed temperature T 
in the sense of thermodynamics and statististical mechanics, even in the presence of the 
gravitational field, when no account is taken of radiation. But we do not know under what 
conditions the temperature will actually be uniform. 

About this question there has been some debate, see e.g. Waldram (1985), page 151. 
It is said that the kinetic energy of each atom in a monatomic gas is 3kT/2 and that, 
when the temperature is the same everywhere, this is paradoxical because it does not 
take account of the potential energy of the atom in the gravitational field. The incident 
involving Loschmidt, Maxwell and Boltzmann has already been mentioned. 



II. 3. The first law and the internal energy 

Is further generalization possible? Can we extend the model to the case that the 
temperature varies with time? The action must be modified, for the temperature becomes 
a dynamical field. Is the temperature one of the variables with respect to which the action 
must be minimized? We need an equation of motion to predict its evolution. (The heat 
equation will be discussed later.) The usual approach is to lay down the additional equation 
by fiat (Section 1.5); is this completely satisfactory? Would it perhaps be preferable to 
have it appear as the result of minimizing the action with respect to variations of the 
temperature field? 

To prepare for the generalization we shall examine some of the main tenets of ther- 
modynamics in the context of the action principle. The question of whether or not it 
is profitable to treat the temperature as a dynamical field variable in the context of the 
action principle can best be assessed a little later (Section III. 3). 

Suppose that the system is in thermal and mechanical isolation except for a force 
that is applied to the boundary. The system is in an equilibrium state with temperature 
T. The applied force is needed to hold the gas within the boundary of the domain E, 
then decreased by a very small amount leading to a displacement of the boundary and an 
increase of the volume by a small amount dV. Assume that this process is reversible. The 
work done by the applied force is 

dW = -pdV. (2.4) 

The first law states that, if the system is in thermal isolation, then this quantity is the 
differential of a function U (T, V) that is referred to as the internal energy of the system. 

Consider the system that consists of an ideal gas confined to a volume V and expe- 
riencing no external forces, not even gravitation. If the gas expands at constant pressure 
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the work done by the gas is pdV and Eq.(1.12) tells us that, 

dV 

pdV = 1ZTM-. (2.5) 

The idea of energy conservation suggests a concept of "internal energy" . It is assumed 
that, under certain circumstances, the work done by the gas is at the expence of an internal 
energy U so that 

pdV + dU = 0, 

or 

UTMdV/V + dU = 0. 

It is an experimental fact (Gay-Lussac 1827, Joule 1850) that the internal energy of an 
ideal gas is independent of the volume (see below) and the more precise statement that 
the internal energy density u is proportional to IZTp is often included in the definition of 
the ideal gas (Finkelstein 1969, page 7). Thus 

u = CyTlTp, U = c v KTM. 

Statistical mechanics gives cy = n, where n is the adiabatic index and takes the value 3/2 
for a monatomic gas. Thus TZTMdV/V + dU = TLTMdV/V + nlZMdT = 0, which implies 
that 

dT = --^-dV, TocV~ 1/n . (2.6) 
n V 

This relates the temperature to the volume and replaces the statement that U is indepen- 
dent of the volume. The calculation from (2.4) onward was done with the understanding 
that M = pV is fixed. 

As we see it, the expression for the internal energy in terms of V and T appears to be 
somewhat ad hoc, derived from external considerations. At the deepest level the concept 
of energy derives its importance from the fact that it is conserved with the passage of 
time, by virtue of the dynamics. The defining equations of hydrodynamics do not admit 
a first integral, and no unique concept of energy; this is a difficulty that our limitation 
to laminar flow, and the action principle, will allow us to overcome. In modern versions 
of thermodynamics, and especially in the thermodynamics of irreversible processes and in 
radiation thermodynamics, conservation laws are all important, but they are postulated, 
one by one, not derived from basic axioms as is the case in other branches of physics, and 
they have a purely formal aspect since they serve only to define various fluxes. See e.g. 
Stanyukovich (1960), Castor (2004). 
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II. 4. The first law and the hamiltonian 



Having adopted an action principle approach we are bound to associate the internal 
energy with the hamiltonian, but one cannot escape the fact that the hamiltonian density 
is defined only up to the addition of a constant multiple of the density. When we decide to 
adopt a particular expression to be used as internal energy over a range of temperatures, 
we are introducing a new assumption. Any expression for the internal energy, together 
with the implication that applied forces increase it by an amount determined by the work 
done, is a statement about a family of systems, indexed by the temperature. This cannot 
come out of the gas law and implies an independent axiom. 

If we adopt the simplest expression for the hamiltonian, 



to serve as "internal energy" , and repeat the analysis of the effect of adiabatically changing 
the volume by means of an applied pressure, then we shall get 



The second term compensates for pdV and so dT = 0, the temperature does not change. 
This is perfectly consistent with the theory as it has been developed so far, but it contra- 
dicts experimental results for an ideal gas. Besides, variation of our present, isothermal 
lagrangian with respect to T does not give a reasonable result, the lagrangian needs to be 
improved. 

II. 5. The adiabatic lagrangian 

In this section we shall treat the adiabatic gas; that is, the polytropic atmosphere with 
Y = 7. But at the end we shall make the case for extending the results to all polytropes. 

The two relations p = IZTp and p = dp 1 imply the relation 1ZT = dp 1 l n between 
the two independent variables T and p that holds for a set of configurations related by 
adiabatic transformations. The index n may be fixed for all configurations, while the 
coefficient d parameterizes the family of adiabats (polytropes). 

It is possible to derive both relations from a principle of least action, by independent 
variation of both temperature and density, but a lagrangian functional of $, p and T can 
only pertain to a single adiabat. An action principle that describes the whole family of 
adiabats must involve additional variables, variables that from a restricted point of view 
of the gas appear as sources of heat. A lagrangian without sources can be interpreted as 
applying to a single adiabat of an isolated system. 

The application that this study is aimed at is the earthly or a stellar atmosphere. 
Those systems are not isolated, but they may nevertheless be treated as being formally 




pdV + dH(T, V) = 0, p = 1ZTM/V. 



In the static case H = KTM\og(M/V) and 



dH = 1ZM log(M/V)dT - TLTMdV/V. 
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isolated when the energy of radiation is included in the hamiltonian and the effect of 
radiation is taken into account through the imposition of boundary conditions. It is felt 
that, if external energy sources are going to be invoked, then it is important that we first 
establish that there is a need to do so. This is why we continue to regard the lagrangian 
as applying to an isolated system. 

Two kinds of additions can be made to the lagrangian without spoiling the equations 
of motion that are essential to hydrodynamics. 

Adding a term independent of p and a term linear in p we consider 

£[<&, p, T] = p($ - v 2 /2 - (p + X) - nTplog(p/po) + pp[T] + f[T\. (2.7) 

The additions do not spoil the relation p = IZTp, nor the continuity of the current. Vari- 
ation with respect to T gives 

pp'[T] - nplog(p/p ) = -f'[T] = -(4a/3)T 3 . (2.8) 

We have set f[T] = (a/3)T 4 , in anticipation of the interpretation of this term as the 
pressure of the photon gas. If the constant a is the Stefan-Boltzmann constant, a = 
7.64 x 10 -15 ergs / K 4 , as it will be taken to be, then this term is very small in most 
circumstances and we must have, on shell, p'[T] 1Z\og(p/ po). The following expression 
will be used 

^[T]=n^Tlog^. (2.9) 
-to 

Eq.(2.8) takes the form 

T n p \ „ , 4a 

-0 P 

and in the important case when n = 3 



^ + log^>+yT3 = 0, (2.10) 



n [ 3 + log lfp-)^ + Y 



which is equivalent to Poisson's law T 3 /p = constant. This reflects the strong affinity that 
is found between the polytropic ideal gas with n — 3 and radiation. The value n = 3 has 
a cosmological significance as well, it is characteristic of the changes in p, p, T induced by 
uniform expansion (Ritter, Emden, see Chandrasekhar 1939, page 48). For other values of 
n, Eq.(2.10) is a mild modification of the polytropic equation of change in the presence of 
radiation. 

The equation of motion that is obtained by variation with respect to p is 

$ - v 2 /2 - (p + X + p[T] = 1ZT(1 + log(p/p ))- (2.11) 
Combined with Eq.(2.10) it reduces, in the static case, to 

p{(j) - A) + (l + log k/k o )KTp = 0, k:=p/T n , (2.12) 
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which has the same form as the equation (1.18) studied in Section 1.5. 

We have thus found an action that, varied with respect to p, $ and T, reproduces all 
the equations that define the ideal, polytropic gas with n = 3. More generally, for any value 
of n, it describes a gas that has its effective polytropic index increased from the 'natural' 
adiabatic value so that the neglect of convection is justified. In the limit when radiation 
becomes negligible it fails to account for the onset of convection, so it should not be used in 
that case. But with this exception it is a mathematical model that incorporates all of the 
physical insight that was summarized at the end of Section 1. 6. 

We suggest that using the lagrangian (2.7), with (2.8) and (2.9), is preferable to the 
usual assumption that (3 := p gas /ptot is constant. Additional evidence for the aptness of 
the adiabatic lagrangian is found in the next subsection. 



II. 6. Energy, pressure and entropy 

The hamiltonian density is, in the static case, with the choice (2.8)-(2.9), 

h = pe + pv 2 /2, pe = # + ^Tplog(^)-^T 4 . (2.13) 

y poT n/ 3 

With the aid of Eq.(2.10) we obtain for the hamiltonian, on shell, when = and v = 0, 

H tot = nllMT + aT 4 V, (2.14) 

in full agreement with the familiar expression for the internal energy of an ideal gas with 
adiabatic index n, augmented by the energy density of the radiation field. This may be the 
first time that this expression for the internal energy has been related to the hamiltonian 
of an action principle. 

The pressure was defined alternatively in terms of the potential, or as the on shell 
value of the lagrangian. We prefer to define the total pressure by the requirement that 

Ptot dV + dH tot = 0. (2.15) 

Taking n = 3 and = we have since T 3 V is constant in this case, 

dHtot = 3RMdT + a(T 3 V)dT = [ - TZMT/V - (a/3)T 4 )]rfV 

and thus 

ptot = TZMT/V + ^T 4 . (2.16) 

This result (2.14-16) is very suggestive. It gives the total pressure as the usual pressure 
of an ideal gas with polytropic index n, augmented by a term that begs to be interpreted 
as a pressure due to heat itself, which is natural when heat is interpreted in terms of 
electromagnetic radiation. Its magnitude is one third of the radiative energy density, as 
expected for the photon gas. 
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An energy conservation equation follows in standard fashion from the action principle, 
namely 

j t (pv 2 /2 + pe) + V- (p^) = 0. 

On the trajectory, p$ = pv 2 /2 + h + p and the standard conservation law results (Castor 
2004, page 14), in the case that there is no contribution from heat flow. 
The internal energy density, with the aid of (2.10) can be expressed as 

pe = nTZTp+ ^T 4 . 

If n is replaced by n', then this is not in agreement with the energy density for a perfect 
gas with adiabatic index n: 

pe = nRTp + -T 4 + (n' - n)TZTp. 

The discreapancy (the last term) indicates that, when n' > n (the case of convective 
stability), if a parcel of air is to move down in the atmosphere, to higher densities, then 
heat has to be supplied at a rate that is proportional to the change in temperature. This is 
precisely the usual interpretation of polytropic change. Namely, one considers a change dT 
in temperature that is accompanied by an addition of heat cpdT. This is a displacement 
along the polytrope with index 

, _ Cp + c 

7 ~ Cy + c' 

or n' = n + c/lZ. Hence 

pe = nIZTp + cpT. 

The increase in energy is thus accounted for; it implies that there is an extra supply of 
energy stored in the gas, distinct from the internal energy and distinct from the contribu- 
tion (a/3)T 4 . This tallies very well with the interpretation of the temperature gradient as 
a greenhouse effect. The extra energy is supplied by the outgoing flux. 



II. 7. Virial theorem 

Both (2.14) and (2.16) are usually derived from considerations outside the proper 
domain of thermodynamics. We prefer an axiomatic foundation of thermodynamics that 
is complete in the sense that it does not need other input. As an example let us discuss 
the use of the virial theorem to make certain predictions concerning stability. 

The virial theorem was introduced into the present context by Kelvin. It is based on 
the scaling properties of the hamiltonian of a system of particles. If H = K + V, kinetic 
energy plus potential energy, then the lagrangian is K — V and the equations of motion 
imply that, up to a time derivative, 
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In the case examined by Kelvin the potential is homogeneous of degree -1, so that, 
in the case of periodic motion , when average is taken over a period, V = 2K. Accord- 
ing to Chandrasekhar (1938) (pp. 49-51), who also quotes Poincare, the internal energy is 
the kinetic energy associated with the microscopic motion of the molecules. It is assumed, 
usually without discussion, that the presence of gravitational forces do not affect the in- 
ternal energy, and that the total energy is obtained by simply adding the gravitational 
potential energy to it. In the present approach there is no place for this argument, the 
hamiltonian is the energy and there is only one energy. 

There is; however, a virial theorem associated with a lagrangian of the type (2.7), that 
we abbreviate as 

C = p($ - v 2 /2) - V. 

(The potential V includes the gravitational field.) Variation of $ and of p give the equations 
of motion 

$ = v 2 /2 + (dV/dp), p = -div(pv), 

which implies that 

If the system goes through a cycle then the average of this quantity over the cycle is zero, 

//* 
dx pv 2 /2 > = < J dx p— > . (2.17) 

In the case of (2.10) we obtain, when n = 3, 

< J dx pv 2 /2 > = < J dx (p((f) - A) + ATZTp + yT 4 ) > . (2.18) 
With Eq.(2.13) this simplifies to 

< J dx pv 2 /2 > = < J dx (p{4> - A) + K(l + log k)pT^j > . (2.19) 

This result, like classical virial theorems, applies exclusively to the case of periodic 
motion. 

In the special case v = Eq.(2.19) is a direct consequence of the equations of motion. 
Such relations, that do not depend on the periodicity of the motion, are not true virial 
theorems. 
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II. 8. The centrifuge and the atmosphere 

Kelvin justified the polytropic model of the atmosphere in terms of radiation and 
convection. Eddington discounted the role of convection and relied on a concept of radiative 
equilibrium. To find out what happens in the case of complete insulation we study the 
analogous situation in a centerfuge. 

Consider an ideal gas. By a series of experiments in which gravity does not play a role, 
involving reversible changes in temperature and pressure, it is found that, at equilibrium, 
the laws p/p = 1ZT and p = kT n are satisfied, constants k, n fixed. When supplemented by 
the laws of hydrodynamics, they are found to hold, or at least they are strongly believed to 
hold, in configurations involving flow, over a limited time span, in the absence of external 
forces. In addition it is said that, at equilibrium, the temperature must be unifoem. 
Keeping an open mind, let us refer to this last statement as "the axiom" . We are talking 
about a fixed quantity of gas contained in a vessel, the walls of which present no friction 
and pass no heat. 

Let the walls of the vessel be two vertical, concentric cylinders, and construct a sta- 
tionary solution of the equations of motion. And why not? We have experimental confir- 
mation of the equations of motion, we applied them to the theory of sound with a degree 
of confidence that is so high that the prediction of rapid variations in temperature may 
never have been subjected to verification (?). In terms of cylindrical coordinates, take 
v z = v r = 0, vg = u, constant. The continuity equation is satisfied with p any function of 
r alone. Then neither T nor p is constant, for the hydro dynamical equations demand that 

ruj 2 = cT\ c = (n + l)1Z ~ 10 7 cm/ sec 2 K (for air). 

At first sight, this seems to violate the axiom, but perhaps not, for this is not a static 
configuration. To save the axiom let us suppose that, by conduction, convection or radia- 
tion, the temperature will tend towards uniformity. Perhaps after a suitably long time has 
passed, T has become constant, in violation of the equations of motion. Let us remember 
that no heat or any other influence is supposed to go by the walls; then surely energy and 
angular momentum must both be preserved during the time that the temperature is level- 
ing out. It also seems reasonable to assume that the final configuration is (macroscopically) 
stationary and uniform, since the existence of fluctuations would imply that the entropy 
had not reached its maximum. But a stationary state with non zero density gradient and 
uniform temperature would seem to contradict the assumptions that we made about the 
gas, which makes the existence of such a state somewhat problematic. 

If we also accept the equivalence principle, then from the point of view of a local 
observer at rest in the flow there is a centrifugal force field, a density gradient and, by 
the laws of Poisson, a temperature gradient. The equivalence principle only applies to 
conditions at one point, and one can question whether the gradient of the temperature 
or of the density is sufficiently local to be covered by the principle. The entire theory of 
relativistic thermodynamics has been founded on the belief that it is (Tolman 1934). 

If we do accept the equivalence principle (without necessarily embracing the tenets of 
traditional relativistic thermodynamics), then we shall be lead to expect that a vertical 
column of an ideal gas, in mechanical equilibrium under the influence of terrestrial gravity, 
and perfectly isolated , will have a pressure and temperature gradient exactly of the form 
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predicted by Lane. This seems to contradict what we think is the prevailing opinion of 
atmospheric scientists, that the temperature gradient owes its existence to the heating 
associated with solar radiation. 

Further measurements in the atmosphere are unlikely to throw light on this, since 
isolation is out of the question. Experiments with a centrifuge may be more realistic. The 
temperature lapse rate is ru 2 x 10 -7 K /cm. If the acceleration is 1000 g at the outer wall, 
then the lapse rate will be .IK /cm. The temperature difference between the inner and 
outer walls will thus be 1 K if the distance is 10cm. In a practical experiment one does 
not have the gas flow between concentric, stationary cylinders. Instead a tube filled with 
the gas is oriented radially on a turntable. Friction against the walls is thus eliminated 
and heat loss is much easier to control. 

On purely theoretical grounds we have come to doubt that complete equilibrium im- 
plies a uniform temperature in all cases. In fact, Tolman (1934, page 314) shows that, 
according to General Relativity, the temperature of an isolated photon gas in a gravita- 
tional field is not quite uniform. The predicted magnitude of this effect is very small, but 
it shows that there are circumstances in which statistical mechanics is not the absolute 
truth. 

The strongest argument that we have found, against the indefinite persistence of a 
temperature gradient in an isolated system, is this. Imagine a large heat bath located in 
the region z > 0. A vertical tube, filled with an ideal gas, has its upper end in thermal 
contact with the bath, otherwise it is isolated. Assume that, at equilibrium, the lower part 
of the tube has a temperature that is higher than that of the bath, then we can run an 
engine by taking out a small amount of heat from the bottom of the tube and returning 
it to the bath. In order for equilibrium to be restored, the heat thus taken out of the 
tube will have to be be restored from the bath, which implies a spontaneous transfer of 
heat in the directio opposed to the gradient of the temperature, in violation of one of the 
statements of the second law (Clausius 1887). According to Graeff (2008), this is exactly 
what he observes! 

If our analysis of the centrifuge is incorrect, what is the mistake? The conventional 
view is that, if a temperature gradient exists at one time, then heat flow will tend to 
dissipate it. Dissipation is determined by Fourier's heat equation; in its general form 

^-T = diC ij djT. (2.20) 

Our error consists of the fact that we are using a lagrangian that contains neither space 
derivatives nor time derivatives of the temperature. We are neglecting phenomena that 
have a long time scale, such as conduction and dissipation . This can be justified if we limit 
our attention to stationary configurations, when T = 0, if the temperature distribution 
makes the right hand side of the equation equal is zero. But a more complete understanding 
of the phenomena requires that we consider nonstationary situations as well, and for that 
we shall have to include additional terms in the action, terms that become significant when 
radiation is not. 

It is not just the paradox of the centrifuge that forces us in this direction. It is 
important that the theory should account for a situation in which the radiation is subject 
to temporal variations, as is certainly the case for the earthly atmosphere. 
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III. Sources 

III.l. Generic source 

We have proposed to extremize the lagrangian with respect to all three fields, $, p and 
T. We have found an expression for the adiabatic lagrangian density, the Euler-Lagrange 
equations of which pass the two tests: 1. When the effect of radiation pressure is neglected 
they give precisely the equations that govern the polytropic atmosphere, for all n. 2. In 
the case that n = 3 and with the radiation term included, they reduce to the equations of 
radiative equilibrium. 

The adiabatic lagrangian describes a single adiabat. To remove this limitation in a 
formal way, let us add another term to the lagrangian density, 

C = p(<i> - v 2 /2 -gz + X)- TZTplog A + f[T] + pTS, (3.1) 

Kq 

where S is an external source. The factor p in the source term is natural and the factor T 
is chosen to make S play the role of a local adiabatic parameter. We have introduced the 
variable k and the parameter /c by 

p = kT n , /9 = W, (fco-1); 

Then k~ 1 l n is Emden's "polytropic temperature". It will be recalled that k parameterizes 
a family of adiabats; in fact, for an isothermal expansion, the variation of — IZlogk is 
precisely the change in specific entropy. The introduction of the source S turns —1Z log ko 
into a field with the interpretation of entropy. We no longer need the parameter and so, 
following Lane, we shall use units of density such that ko = 1. 

The internal specific entropy is 1Z log(T n /p) and the total specific entropy is 

Stot = ft log — + S. 

p 



With this convention 



C = p(<& - v 2 /2 - 4> + A) + pTS tot + ^T 4 , (3.1) 



where </> is the gravitational potential. Variation with respect to T leads to 

r) An 
P -^{TS tot ) + -T* = Q. (3.2) 

As an equation for S to t it has the general solution 

Stot = -^-\v[p\. (3.3) 
3p pT 

Taking this as the definition of the potential (V is the value of pTStot + f T A at the 
extremum with respect to variation of T) we have 

C = p(<5>-v 2 /2-<p + \) -V[p\. (3.4) 
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The gradient of the equation obtained by variation of p is 

~P-jX°- pgrad^ = gradp, (3.5) 

with 

P = P V'-V=(1- P^ p )(pTS tot + t|t 4 ) = ~Tp 2 ^ + a -T\ (3.6) 

The last equation is justified by the fact that the partial derivative of pTS to t + (a/3)T 4 
with repect to T vanishes, Eq.(3.2). 

We shall verify some important relations of thermodynamics, and for this we must 
take T and p to be constant, with M = pV, and = 0. In this case 

p = MT^^ + -T 4 = UMT/V+-T 4 + MT^-. (3.7) 
oV 3 3 oV 

The hamiltonian density is, in the static case, in the absence of gravity, 

h = -pTStot ~ implying that U = -MTS tot - ^T 4 V. (3.8) 
Variation of h with respect to T gives zero on shell, so this is the same as 



Thus 



u = (1 - Td T )h = pT 1 -^- + aT = nKpT + aT 4 + pT 2 — 

U = MT 2 -^+aT 4 V = n1lMT + aT 4 V + MT—. (3.9) 
Using (3.7) and (3.8) one verifies that 

^ = (Td T - l)p, (3.10) 

an important consequence of the existence of entropy in general. See Finkelstein (1969) 
page 26. Also, 

dU = MT-^-dT + Tdi^MT-^j + 4aT 3 VdT + aT 4 dV, 

pdV = MT^^dV + -T 4 dV, 

and the sum is dU + pdV = dQ = T(dp/dT) = 0, the last on shell. 

We used the last expressions in (3.7) and (3.9) because they are familiar, but if we 
return to (3.8) and the first expression for p in (3.7) we see immediately that dU+pdV = 0. 

If instead we consider a change that involves outside forces acting via the source, then 
dU+pdV — MT5S, which confirms the interpretation of S as a contribution to the specific 
entropy. 
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The calculations that have been presented in this subsection are offered as proof 
that the variational approach that is being advocated is fully compatible with classical 
thermodynamics. This gives us faith in the basic framework and courage to proceed on a 
more speculative course. 



III. 2. Electromagnetic fields 

We write the Maxwell lagrangian as follows, 

£rad = —D 2 - -H 2 + D ■ (dA -A)-H-dAA + JA, (3.13) 

and add it to the ideal gas lagrangian 

£ gas = p($ - v 2 /2 -cf> + \)-KTp log k + ^T 4 , (3.14) 

Since the susceptibility of an ideal gas is small, the dielectric constant may be expressed 
by 

e = l + K[p,T], or - = l- K [p,T}. (3.15) 

Paramagnetic effects will be ignored at present. An interaction between the two systems 
occurs through the dependence of the susceptibility on p. The source S has become 
— (D 2 /2p)(k/T). If this quantity has a constant value then it produces a shift in the value 
of the parameter k. 

Two interpretations are possible. The electromagnetic field may represent an external 
field, produced mainly by the source J, and affecting the gas by way of the coupling implied 
by the dependence of the dielectric constant on p. Alternatively, J = and the field is 
produced by microscopic fluctuations, quantum vacuum fluctuations as well as effects of the 
intrinsic dipoles of the molecules of the gas. In this latter case the main effect of radiation 
is represented by the radiation term aT 4 /3. Our difficulty is that neither interpretation 
is complete, and that we do not have a sufficient grasp of the general case when either 
interpretation is only half right. The following should therefore be regarded as tentative. 

Variation of the total action, with lagrangian £ r ad + C gas , with respect to A, D, H and 
T gives 

D = dAH, (3.14) 

A = D/e, (3.15) 
/iH = -dAA, (3.16) 

and 

n(n-\ogk)p-\^ + A ^ = Q. (3.17) 
Taking into account the first 3 equations we find for the static hamiltonian 
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With the help of (3.17) it becomes 



H = J d*x(cj>p + nllpT +^ + ^ + aT 4 ) - J d*x T^-^. (3.18) 

The last term, from the point of view of the thermodynamical interpretation of electro- 
statics, is recognized as the entropy (Panofsky and Phillips 1955). On a suitable choice of 
the functional k it merges into the internal energy. For example, if k = pT it takes the 
form pTS with S = D 2 . 



III. 3. Discussion 1. Using T as a dynamical variable 

The idea of extremizing thermodynamical potentials with respect to the temperature 
is far from new, but in the context of the action principle it is likely to raise questions. 

Radiation Hydrodynamics (Castor 2004) is defined by the following equations (without 
specialization to laminar flow) : the continuity equation, the hydrodynamical equation, the 
equation of motion, and an 'energy equation'. This last equation takes the place of the 
equation that results from variation of the temperature. To see this we have only to review 
the canonical conservation of energy in lagrangian/hamiltonian form. 

With C as in (3.4), 

C = p($ - v 2 /2-cf> + A) -V[p,T], 
with no derivatives in the functional V[p,T]. We have 

dC .dC -dC dC -dC 

— = p h $ — - + grad$ ■ h T — . 

dt 'dp d§ <9grad$ dT 

On shell, the first term on the right is zero. The equation can be rearranged to read 

Dh . . -dC 

_ + dlvM = -T-. 

The 'energy equation' is the vanishing of the left side; it is thus equivalent to setting 
dC/dT = 0. There is no evidence of a contribution to energy transport by heat flow, but 
that is not a problem if we accept Emden's postulate of a constant heat flow. 

The effect of heat transfer by radiation and by conduction has been left out, for 
simplicity and since it does not immediately affect the polytropic atmosphere. But it 
means that our treatment is incomplete. Most important, it prevents us from dealing with 
the isothermal atmosphere that is, after all, an important special case. Consequently, it is 
felt that a general understanding of the effect of radiation is lacking. 

Let us examine the total lagrangian, 

£ = £ rad + £ gas = p($ - v 2 /2 - cj) + A) - TZTp log ^ 
+ — + ^H 2 + D ■ (8A -A)-H-dAA + JA+ -T\ (3.19) 
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So long as e, p, and J are independent of p, T and v, the variational equations of motion 
that are obtained by variation of v, p, A, H and D are all conventional, at least when n = 3 
(for all n if radiation is neglegible). It would be possible to be content with that and fix T 
by fiat, as is usual; in the case of the ideal gas without radiation the result is the same. But 
if e depends on p and on T, which is actually the case, then we get into a situation that 
provides the strongest justification yet for preferring an action principle formulation with T 
as a dynamical variable. The equations of motion include a contribution from the variation 
of e with respect to p, so that one of the basic hydrodynamical equations is modified. Thus 
it is clear that the extension of the theory, to include the effect of radiation, is not just a 
matter of including additional equations for the new degrees of freedom. The presence of 
the term D 2 /2e[p, T] certainly introduces the density p into Maxwell's equations; that it 
introduces D into the hydrodynamical equations is clear as well. The over all consistency 
of the total system of equations can be ensured by heeding Onsager's principle of balance, 
but the action principle makes it automatic. 

Variation of the action with respect to T offers additional advantages. The usual 
procedure, that amounts to fixing p = kT n , k and n constant, gives the same result when 
radiation is a relatively unimportant companion to the ideal gas, but in the other limiting 
case, when the density is very dilute and the gas becomes an insignificant addition to the 
photon gas, it is no longer tenable. We need an interpolation between the two extreme 
cases and this is provided naturally by the postulate that the action is stationary with 
respect to variations of the temperature field. 

In the absence of the ideal gas we have another interesting system, the pure photon 
gas. The analogy between the photon gas and the ideal gas is often stressed; there is an 
analogue of the polytropic relation that fixes the temperature in terms of p; the pressure 
of the photon field is (a/3)T 4 . Our lagrangian already contains this pressure; we should 
like to discover a closer connection between it and the electromagnetic field. In the limit 
when the density of the ideal gas is zero, Eq.(3.17) becomes 

D 2 Ok 4a , 

H T = 0. 



2 dT 3 

In the absence of the gas it is reasonable to impose Lorentz invariance, so we include 
magnetic effects by completing the last to 



2 dT 3 

If we suppose that k[p, T], in the limit p = 0, takes the form <xT 2 , then 

aF 2 = — T . 
3 

The radiation from a gas of Hertzian dipoles can be shown, with the help of the Stefan- 
Boltzman law and Wien's displacement law, to satisfy a relation of precisely this form. 
Whether the same relation holds in vacuum is uncertain, but it is suggested by an analysis 
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of the effective Born-Infeld lagrangian calculated on the basis of the scattering of light by- 
light (Euler 1936, Karplus and Neuman 1950). See also McKenna and Platzman (1962). 

III. 4. Discussion 2. The temperature gradient of the atmosphere 

We return to the question of the heating of the atmosphere by solar readiation. The 
susceptibility of air may be approximated by k = ap/T (a constant); in which case the 
lagrangian includes the following 

-TZTp log p + nRpT log T + apT~ 1 D 2 /2. (3.20) 

Up to a factor pT, the first term is the entropy of the isothermal gas. The first two terms 
together give the entropy of the polytropic gas. The entropy introduced by radiation is 
aD 2 (8(kT) / dT) , namely zero under the assumptions made, which would suggest that the 
greater effect of radiation on the specific entropy of air comes about through the influence 
of the term aT 4 /3 in the lagrangian. The second term was included in our theory in 
order that variation of the action give both of the Poisson relations, the ideal gas law 
and p oc T n . Let us return to the problem of determining the origin of the temperature 
gradient, the physical reason for the fact that the atmosphere is nearly adiabatic. In 
the model, adiabacity comes with the inclusion of the term riRpT log T in the action. If 
this term could be linked with the last term in (3.20); that is, if it could be shown that 
aD 2 /2 m nTZT 2 logT, then we could conclude that the temperature profile is created by 
the radiation. But this seems unlikely, as it amounts to fine tuning. If we suppose that 
the term in question is responsible for the bulk of the effect of solar radiation on the 
atmosphere of the earth, then we must admit that it undergoes important diurnal and 
seasonal variation. In fact, the polytropic index of the troposphere is usually quoted as 
a constant, without any indication that important variations have been observed over a 
period of time. 

It could be argued that the relatively stable polytropic index of the earth is evidence 
that the atmosphere, if isolated, would continue to manifest a temperature gradient and 
that an ideal gas, isolated in a gravitational field, may not tend to an equilibrium state with 
uniform temperature, but the explanation of the temperature gradient in terms of convec- 
tion is more attractive by far. See the last paragraph of Section 1.6. Additional terms 
(containing in particular derivatives of the temperature) are needed in the lagrangian, al- 
though they are insignificant in applications to the polytropic atmosphere under conditions 
of strong radiation fields. 



27 



IV. Stability of atmospheres 

IV. 1. The isothermal column 

We consider the space that is tangent to a static solution with density pq. Setting 
p = Po + S p we have the following equations for the perturbation Sp, 

-v = KT(8p/p o y, Sp = -(p v)', 

where the prime denotes differentiation with respect to z. Thus 

6p = lZT(p ay, a:=Sp/p . (4.1) 

For a harmonic mode with frequency u, —u 2 Sp = RT(poa')' , and 



u 2 



TIT 



J a 2 podz = J a(poa')'dz = 5pa' —J p a' 2 dz. 



The configuration is stable if this implies that uj 2 > 0, which will be the case if the 
boundary term vanishes. To justify any choice of boundary conditions we have only the 
conservation of mass, f Spdz = 0. This ensures that Sp fall off at infinity and we are left 
with -6p(0)a'{0). 

We shall show that a'(0) = 0. Eq.(4.1) tells us that 

uj 2 

;« = (Po/PoK + a" = -^a' + a". 



TIT iru,ru; TIT 
This is a linear differential equation with constant coefficients, with general solution 



a = Ae k + z + Be k - z , k ± = ± W (^-) 2 - — 



2TZT V x 2nT' m 
Since, up to an irrelevant constant factor, p = exp(—gz/7ZT), 



Sp = Po a = Ae a+Z + Be a - Z , a± = — J= ± J i^;) 2 - £=. 

These functions are integrable only if < uj 2 < g 2 /41ZT, and in that case 

f e i A B A B 

5M = / Spdz = = — + — , 

J a+ a_ k- K+ 

the vanishing of which requires that a'(0) = 0. When u 2 > g 2 /ATZT we have instead 
to do with a contour integral, and reach the same conclusion. Therefore, not only is the 
condition u 2 > verified; it is also confirmed that the boundary condition a'(0) = is the 
only one possible. We have seen that this choice of boundary conditions is the one that 
ensures the conservation of mass. 
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IV. 2. The polytropic column 

Let us leave the parameter k = p/T n free and fix the value of n. This conforms to 
the usual approach when the temperature is fixed by edict, but it is consistent with our 
formulation if n = 3 only. We study the stability to vertical perturbations. 

The static solution is 

cT = X-gz, c:=7£(l + logfc). 
A first order perturbation satisfies 

<5> + 5X = c5T, thus v = -c5T', (JV.l) 

and 

p=-(vp)', p = -(v p )\ and nT^dT = c(T n 5T')'. (IV.2) 

Let x = X/g — z, < x < X/g and let /' = df jdx from now on. Solutions of the type 
ST = exp(iut) f (x) satisfy the equation 



(x n ST')' + —(x n 5T) = 0, v 2 = nu 2 /g. (JV.3) 



, 2 

X 

The solution that is regular at the origin of x (the top of the atmosphere) is 

5T = oF 2 (n,-v 2 x)e iujt . 

The generalized hypergeometric function is positive for positive argument and it oscillates 
around zero for negative argument. 

Boundary conditions . If we fix ST = at the bottom of the column we can prove stability 
as follows. For a harmonic perturbation, 

v 2 J x^iSTfdx = - J (x n 5T')'5Tdx = J x n {5T') 2 dx > 0, 

which shows that v 2 is positive and that the solutions are oscillatory in time. But there is 
no justification for this choice of boundary condition. 

It is not unusual to fix the upper boundary, and to require that the perturbation 
vanish there. If v 2 is positive the argument of the hypergeometric function is negative. 
The function oscillates around zero and for a discrete set of values of the frequency it 
vanishes at the upper end. If v 2 is negative then the hypergeometric function is positive, 
and the boundary condition cannot be met. 

However, there seems to be no better reason to fix the upper boundary. The natural 
boundary condition is that the mass must be preserved, thus 

5M = J 5pdx = J T n - x 5Tdx = 0. 
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This may happen for a discrete set of positive values of v 2 . For negative values of v 2 the 
integrand is definite so that it can not happen. The calculation is valid only in the case 
n = 3; this atmosphere is stable. For other values of n the calculation is more difficult. 

The problem can be converted to a standard boundary value problem by rescaling of 
the coordinate. 

The mass is 



M = Ak(-y dxx 6 = — (-Y-. (IVA) 
a J 4 a g 

The proper definition of gravitational energy is ambiguous, but (11.21) suggests that it is 

E g = Ag [ p(gz-\)dz = -Akg( 9 -f ^ dxx' = -^(-f- 
J a Jo 5 a g 

The last expression, and those that follow, refer to the static solutions. The thermodynamic 
part of the hamiltonian is 

r 2 dr(lZTp log k + V) = / rfxT 4 
3 4 J R 

The integrand on the right hand side of (11.21) is thus E g + 4(H — E g ) = 0, as it must be. 



IV. 3. The polytropic gas sphere. The hamiltonian 

Here we study the self gravitating polytropic gas. A correction is needed in the 
expression for the lagrangian, and we need to take care with respect to the definition of 
the gravitational potential. 

First of all, it would not be difficult to argue that the correct expression for the 
gravitational energy is 

-% [ d^x' P ^'\ (IV.6) 
2 J \x x \ 



Eg 



and that, consequently, the term —gz in the lagrangian has to be replaced by —E g . Nothing 
else is needed, but to make an important point it will be useful to introduce a potential, 
a functional cj)[p] defined by 

<P[p](x) = 0(f) = / d 3 x^^± + 0[O]. (IV.7) 
J \x — x I 

The last term, <p[0] := (p[p]\ p= o is of course an arbitrary constant field. The value chosen 
for this constant is irrelevant, but it must be kept in mind that it is chosen once and for 
all and that it is independent of p . The sign is opposite to that used by Eddington; it is 
chosen so that the gravitational force is -grad (p. 
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then 



For any spherically symmetric distribution define 

M(r) =4tt ! r' 2 p{r')- 
Jo 

C r C KA 

0(f)^0(r)= / -L^dr' + m, 
Jo r 

r 2 (j)' = GM, M' = 4nr 2 p. 



(IV.8) 



(IV.9) 
(IV.10) 



Since first offered by Perry (1899) it has become customary to present a criterion for the 
stability of static configurations, based on an evaluation of the energy. The better to 
understand it we replace the definition (4.6) of the gravitational energy by 



E, 



, - \ I P<t>d 3 x 



\ J M'<j>dr, 



(JV.ll) 



acknowledging that, since 4> is defined up to an additive constant, the same is true of the 
energy. 

As we have seen in Section II. 7, it is reasonable to identify the total energy (including 
the gravitational energy and the internal energy) with the hamiltonian. The lagrangian is 



a. 



L = p($ — v 2 /2 - cf>/2 + A) - TZTp log k - -T . 



{IV. 12) 



The factors of 1/2 in the term p<p/2 and in (IV. 11) arise from the fact that p<p is a 
homogeneous functional of p of order 2, 5(p4>/2) = (p. Variation of the action with respect 
to T gives the relation 

An T 4 

KT{n - logfc) = . (JV.13) 

3 p 



The hamiltonian density is 



h = pv 2 /2 + + TZTp log k + ^T 4 , 



or, in view of (IV. 13), 



h= 1 -v 2 + 1 -p ( j ) +^Tp, c = nq + iogk), 

whence the hamiltonian (= total energy) 

H = J d 3 xh = E g + J d 3 x(^v 2 + ^Tp). 

In the static case we have the equation of motion 

A = (p + cT. 
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(IV.U) 



(JV.15) 



Since the surface of the star is at the point where T = 0, it follows that 

A = <p(R) = 0(0) + cT(0). (IV.16) 
Applying J d 3 xp to (IV. 15) we obtain 

MX = 2E g + J cTpd 3 x. 

Hence (4.14) reduces, in the static case, to 

H = E g + \(M<f>{R) - 2E g ) = -E g /2 + \m${R). (JV.17) 

Note that H = E g + U and that this is what Eddington and others call the total energy. 

It remains to calculate E g , and here we follow Eddington. To begin, using only the 
definitions (4.8)- (4.11), 

E g = 4tt J ]^p(j)r 2 dr =^J M'<f)dr = ^M(j>(R) - ^ J M<p'dr, (4.18) 

and 

-\ J Mfdx = -i / ^Idr = \ J GM\\)'dr =~~f jMdM. (4.19) 
Next, the poly tropic relation can be used to show that 



^ J M<f>'dr = J jMdM. 



Hence / M(p'dr = 2GM 2 /R and finally 

1 CM 2 1 1 CM 2 

E g = -Mcf>(R)-—, H =-M<KR) + - — . (4.20) 

So far, the only difference between our calculation and those of Eddington and others is the 
fact that we have left open the zero point of the potential. Eddington's field is 0(Eddington) 
= —(p — GM/R and his boundary condition ^(Eddington) = at the surface, amounts to 

4>{R) = — — 55— • (Eddington's choice) 
R 

According to (4.16) this is the same as <p(0) = —(MG/1Z) — cT(0). But, as we have 
emphasized already, <p(0) = <p[0] is a constant functional, independent of the dynamical 
variables. Therefore, Eddington's choice is not only ad hoc but, in the context of the 
action principle, wrong! In fact, we know of no physical theory in which the manifold of 
physical states is restricted to a single energy surface in phase space. 
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The total energy provided by the action principle is given by (IV. 20) and (IV. 16), 



H[p,T}= 1 -M(^f + cT(0) + ct>[0ty 



It depends on the value M chosen for the constant of the motion, the initial value T(0), 
and the choice of an inessential zero point for the gravitational potential. 

We feel justified to conclude that the insistence on an action principle is much more 
than an aesthetic preference; it is an essential aid to avoid fortuitous conclusions. 

We have chosen to investigate the case n = 3, since Eddington's calculations are 
valid in that case only. As was shown, they lead to no conclusion even in that case. The 
statement that the static configurations are stable for n < 3 and unstable for n > 3 may 
be correct, but to say that it was proved by the argument first advanced by Perry, or by 
the same calculation repeated in many of our modern textbooks, is an exaggeration. 



IV. 4. The polytropic gas sphere. Stability 

We use the lagrangian 

£ = p($ -v 2 /2-(j)/2 + \) - nTplogk+-T 4 , k:=p/T n . (JV.21) 

3 

Variation with respect to T gives 

n(logk-n) = ^-T 3 /p. (IV.22) 

With n = 3 this makes k a constant, and log/c = 3 when radiation is neglected. In the 
remainder of this section, we set, for all values of n, 

p = kT n , k constant. 

This is the usual polytropic relation used by Eddington and others, but it is consistent 
with (4.21) only when n = 3. The remaining dynamical equations are 

UV = 0' + cT', p + r- 2 (r 2 pv)' = 0, 



Dt 

4nGp = r- 2 (r 2 (P'y, p = kT n . 

• The static solution. Eliminate <j) by 4>' = —cT' and change variables, setting r = x/a, 
a constant, Poisson's equation becomes 

^^x 2 T n + (x 2 T')' = 0, 
ca 2 
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where the prime now stands for differentiation with repect to x. Set f(x) = T(x)/T(0) 
and a = ^4nG/cT(0) so that finallly 

x 2 f n +(x 2 f')' = 0, /(0) = 1, /'(0) = 0. 

The solution decreases monotoneously to zero at x = X, this point taken to be the surface 
of the star. At the outer limit f(x) oc X/x — 1 + o(X — x) n . The integration is done easily 
and accurately by Mathematica, especially so for integer values of n. The radii are, for 
n = 2 : X = 4.355, n = 3 : X = 6.89685636197, n = 4 : X = 14.9715. 

• For the fluctuations we assume harmonic time dependence, then the equations are 



-u 2 r 2 6p = (r 2 p (8(f)' + c5T'))', bp = nkT^ST, 

4nGr 2 5p = (r 2 5(P'y. 
Introduce the function 5A4 = r 2 5(j)'. Eq.s(IV.23-4) then take the form 



(JV.23) 
(JV.24) 



4tvG 



5M = p5M + r 2 pcST' + constant, 



where the constant can only be zero, and 

(47rG)r 2 (?iA;T n - 1 5T) = 6M\ 

Elimination of 5T leads to 

SM' 



4nG 



SM = pSM + 



AnGkn 



r 2 p 



x 2 T n ~ 



Changing the scale as before we get 



4nGkT n (0)' 



(JV.25) 



The crucial point is the choice of the correct boundary conditions, at x = as well as 
the outer surface (x = X). At the center the solutions take one of two forms, 1 + Cx 2 + 
which is unphysical, or else x 3 + Cx 5 + ... Accordingly we set 



5M(x)=x 3 g(x), (/(O) = 1, g'(0)=0. 



(JV.26) 



The boundary conditions at the outer boundary are determined by the fact that the mass 
is conserved, 

5M = 5M(X) = 0. 

The equations then imply that the null point is of order n. With these boundary conditions 
(IV. 25) becomes a well defined Sturm-Liouville problem with an essentially self adjoint, 
second order differential operator. 
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Numerical calculations with the help of Mathematica are not difficult in the case of 
integer values of n. It is found that, when n = 2 and for n = 3, 5Ai(X) is positive in 
the whole range, for all negative values of v 2 and for positive values below a limit v 2 that 
is about .06 for n = 2 and compatible with for n = 3. The latter is the first, nodeless 
solution of a sequence of solutions that we have not determined in detail. The function 
falls to zero at the surface, where there is an nth order zero. Above this lowest value of v 2 
is a discrete set of other values of v 2 at which the boundary condition is satisfied. 

At the special value n = 3 the 'ground state', the lowest value of z/ 2 , has approached 
very close to zero. 

Polytropes with n = 4 are widely believed to be unstable, but a positive proof of 
this is not known. We have searched for harmonic solutions with negative values of v 2 . 
The value n = 4 is indicated because it is the only integer in the interesting range, and 
because Mathematica is much more managable in this case. (Accuracy is lost when non 
integral powers of negative numbers appear at the end point.) There seems to be a discrete, 
decaying nodeless mode with v 2 = —.015796, but a bifurcation at this point in parameter 
space makes the conclusion uncertain. We carried the calculation to 15 significant figures 
in v 2 but solutions do not converge towards a function that vanishes at the surface. To 
overcome this difficulty we reformulated the problem in terms of the variational calculus. 
The "solution" found for v 2 = —.015796, truncated near both ends, was used as a trial 
function, to show conclusively that the spectrum of v 2 extends this far. Among may papers 
on this topic we mention Cowling (1936) and Ledoux (1941). 



IV. 5. The case n = 3 

This case is widely believed to mark the boundary between stable and unstable poly- 
tropes. The equations are conformally invariant and a time independent solution is found 
by an infinitesimal conformal (homology) transformation, 

df = rf + f. (IV.26) 

This does not represent an instability, but a "flat direction" , a perturbation from which the 
system does not spring back, nor does it run away. There must also be a second solution, 
linear in t, of the form 

Sf = t(rfY, 5p = t(rp' + 3p). 

The equation of continuity becomes rp' + 3p + vp' + r~ 2 (r 2 v)' p = 0, whence v = — r. 

This linear perturbation is the first order approximation to the exact solution found 
by Goldreich and Weber (1980), of the form 

f(r,t) = ^y/O^ x = r/a(t). 

The continuity equation is solved by v = ax; thus $ = — (d/a)(r 2 /2), and 

- v 2 /2 = -aax 2 /2 = cT + (p. 
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This leads to 

4> = a(t)4> oc / + Ka 2 ax 2 /6, K = 3k 1 ^ 3 /c, 
and Poisson's equation becomes 

f 3 + -^(x 2 f'Y = ^a 2 dx 2 /6 = -Ka 2 d = A, constant. (IV.27) 

X 2 X 2 

There is a first integral, 

— a 2 — \/a = C, constant. 
Rescaling of t and a reduces this to one of three cases 

a = y/l + 1/a, a = y/l — 1/a, a = 1/ y/a, 

but only the first is compatible with analyticity at t = 0, thus 

t = y/a^/l + a — arcsinhv^a. 

Setting a = 1 + b we find 

t= ^/T/2(b-b 2 /2) + o(b 3 ) 

The factor a(t) is zero at a finite, negative value of t and increases monotoneously to 
infinity, passing through 1 at t = 0. We can of course reverse the direction of flow of t to 
get collapse in the finite future. 

Eq. (IV.27) was solved numerically (Goldreich and Weber, 1980). The solution is 
similar to the solution of Emden's equation, just prolonged a little at the outer end, so 
long as < A < .00654376. For larger values of A the distribution does not reach zero 
and increases for large r. For simillar studies of collapsing, isothermal spheres see Hunter 

(1977) and references therein. 

It is sure, therefore, that the polytrope with n = 3 is not stable. Suitably erturbed, 
the star may expand or collapse, until the higher or lower density causes a change in the 
equation of state. Among may paers on collapse we may mention Arnett (1977), Cheng 

(1978) , Hunter (1977) and Van Riper (1978). 

V. General Relativity 

V.l. Lorentz invariance 

The limitation to small velocities, small compared to the velocity of light, is justified 
almost always, with the sole exception of the photon gas. We shall now modify our treat- 
ment of the non relativistic gas of massive particles to make it consistent with relativistic 
invariance. 

We need a 4-dimensional velocity and an associated velocity potential, 

= d^ip =: ip^, ^ = 0,1,2,3, 
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where ip is a scalar field. There is only one reasonable lagrangian (Fronsdal 2007), 

C= P -{g^iP,^, v -c 2 )-V[p\. (5.1) 

The metric is the Lorentzian g = diag (c -2 , —1, —1, —1). In the case of velocities small 
compared to c we set 

if; = c 2 t + $ 

and find to order o(c -2 ) the non relativistic lagrangian (1.10). Henceforth c = 1. 
We easily allow for a dynamical gravitational field by generalizing the measure, 

A = J dtd 3 xy^gC. 

In a weak, terrestrial gravitational field the usual approximation for the metric is 
g = diag. (1 — 2gz, —1, —1, —1), which leads to (2.10). 

The concept of energy (density) is all-important in thermodynamics and in relativistic 
field theories but ill defined in General Relativity. However, as long as we limit our atten- 
tion to time independent configurations, we expect to be on relatively safe grounds when 
we identify the energy density with the time-time component of the energy-momentum 
tensor, 

In the non relativistic limit Too is our hamiltonian augmented with the rest mass. 
The Euler-Lagrange equations include the conservation law 

d^ = 0, J" := V^g^u. 

The integral J y/^gpd 3 x is a constant of the motion (for appropriate boundary conditions) 
and can be interpreted as mass. This is an essential improvement over the traditional treat- 
ment. A conserved current also permits an application to a non neutral plasma (Fronsdal 
2007). The (conserved) mass plays a central role in fixing the boundary conditions in the 
non relativistic theory; to retain this feature in the relativistic extension is natural. 

V.2. Polytropic star with radiation 

Here we propose to try out the lagrangian (2.7) or its relativistic version for the 
mixture of an ideal gas with the photon gas. In the case that the radiation pressure is 
relatively unimportant there is nothing new in this, and in the special case that n = 3 the 
theory is identical with that of Eddington. 

In the relativistic case the action principle offers advantages even in this particular 
case. Clarification of the role of mass, which is confused or at least confusing in the 
traditional treatment, is an important part of it. Another advantage is the relative ease 
with which one may proceed to study mixtures. 

Variation of the action with respect to the temperature gives the relation (11.10) that 
shows a departure from the polytropic relation p = kT n when n^3. (If this last relation 
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is accepted, in lieu of (11.10), then from this point on the equations of motion are the same 
as with other methods.) The relation between Eddington's parameter f3 and k, n is 

1 — / - 1 a 

^-Ptot/Pgas-l + ^p 

It is constant only when n = 3. In the relativistic theory, the same relations hold; Eq.(II.lO) 
remains valid. The equation that determines the temperature is transcendental; the first 
approximation is logk/ko = 3. 

Applications to real stars should await the incorporation of heat flow, not important in 
the case of an isolated atmosphere and of secondary importance in the case of the earthly 
atmosphere, but perhaps vital for a full understanding of the physics. 



VI. Conclusions 

VI. 1. On variational principles 

Variational principles have a very high reputation in most branches of physics; they 
even occupy a central position in classical thermodynamics, see for example the authori- 
tative treatment by Callen (1960). An action is available for the study of laminar flows in 
hydrodynamics, see e.g. Fetter and Walecka (1960), though it does not seem to have been 
much used. Without the restriction to laminar flows it remains possible to formulate an 
action principle (Taub 1954, Bardeen 1970, Schutz 1970), but the proliferation of velocity 
potentials is confusing and no application is known to us. Recently, variational principles 
have been invoked in special situations that arise in gravitation. 

In this paper we rely on an action principle formulation of the full set of laws that 
govern an ideal gas, in the presence of gravity and radiation. To keep it simple we have 
restricted our attention to laminar, hydrodynamical flows. 

It was shown that there is an action that incorporates both of Poisson's laws as 
variational equations, the temperature field being treated as any other dynamical variable. 
The idea of varying the action with respect to the temperature is much in the classical 
tradition. The variational equations of motion are exactly the classical relations if radiation 
is neglected, or if n = 3. 

The first encouraging result comes with the realization that the hamiltonian gives the 
correct expression for the internal energy and the pressure, including the contributions of 
radiation, under the circumstances that are considered in classical thermodynamics; that 
is, in equilibrium and in the absence of gravitation. This is an indication that the theory is 
mathematically complete, requiring no additional input from the underlying microscopic 
interpretation. This conclusion is reinforced by an internal derivation of a virial theorem. 

Into this framework the inclusion of a gravitational field is natural. Inevitably, it 
leads to pressure gradients and thus also temperature gradients. If other considerations, 
including the heat equation, are put aside, then the theory, as it stands, predicts the 
persistence of a temperature gradient in an isolated system at equilibrium. The existence 
of a temperature gradient in an isolated thermodynamical system is anathema to tradition, 
and further work is required to find the way to avoid it, or to live with it. Physical 
considerations indicate that the answer is to be found in the phenomenon of convection. 
The theory in the present form can be applied when convection is not important. 
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A secondary but satisfying result of this work has been the application of the action 
principle to the study of the energy concept. Without a well defined hamiltonian it is quite 
impossible to attach an operative meaning to any expression for the value of the energy; 
it is always defined up to an additive constant, independently for each solution of the 
equations of motion. With a hamiltonian at our disposal we are in a position to give voice 
to our misgivings concerning the way that "energy" has been invoked in some branches of 
physics over a period of over 100 years. Though we conclude that past demonstrations of 
instabilities of polytropes are inconclusive, we do not suggest that the results are wrong. 
It is of course agreed that n = 3 represents an important bifurcation point. 

We have insisted on the role played by the mass in fixing the boundary conditions, ver- 
ified for 3 different atmospheres. The existence of a conserved current and the associated 
constant of the motion is especially important in the context of General Relativity where 
the absence of this concept casts a shadow of doubt on the choice of boundary conditions 
(Fronsdal 2008). Indeed it is strange that the equation of continuity, a major pillar of non- 
relativistic hydrodynamics, has been abandoned without protest in the popular relativistic 
extension. See Kippenhahn and Weigert (1990), pages 12-13. 

The interaction of the ideal gas with electromagnetic fields has been discussed in a 
provisional manner. The transfer of entropy between the two gases is in accord with the 
usual treatment of each system separately. 

VI. 2. Suggestions 

(1) It is suggested that observation of the diurnal and seasonal variations of the equa- 
tion of state of the troposphere may lead to a better understanding of the role of radiation 
in our atmosphere. The centrifuge may also be a practical source of enlightenment. We 
understand that modern centrifuges are capable of producing accelerations of up to W 6 g. 
Any positive result for the temperature gradient in an isolated gas would certainly have 
important theoretical consequences. 

(2) We suggest the use of the lagrangian (2.7), or its relativistic extension, with T 
treated as an independent dynamical variable and n' = n. Variation with respect to T 
yields the adiabatic relations between p and T, so long as the pressure of radiation is 
negligible, but for higher temperatures, when radiation becomes important, the effect is to 
increase the effective value of n' towards the ultimate limit 3, regardless of the adiabatic 
index n of the gas. See in this connection the discussion by Cox and Giuli (1968), page 271. 
In the case that n = 3 there is Eddington's treatment of the mixture of an ideal gas with the 
photon gas. But most gas spheres have a polytropic index somewhat less than 3 and in this 
case the ratio (3 = p gas /ptot m ay not be constant throughout the star. The lagrangian (2.7), 
with n identified with the adiabatic index of the gas, gives all the equations that are used to 
describe atmospheres, so long as radiation is insignificant. With greater radiative pressure 
the polytropic index of the atmosphere is affected. It is not quite constant, but nearly 
so, and it approaches the upper limit 3 when the radiation pressure becomes dominant. 
Eddington's treatment was indicated because he used Tolman's approach to relativistic 
thermodynamics, where there is room for only one density and only one pressure. Of 
course, all kinds of mixtures have been studied, but the equations that govern them do not 
supplement Tolman's gravitational concepts in a satisfactory manner, in our opinion. Be 
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that as it may, it is patent that the approximation (3 = constant, in the works of Eddington 
and Chandrasekhar, is a device designed to avoid dealing with two independent gases. 
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